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ABSTRACT 

We report on the discovery and analysis of the transient X-ray pulsar CXOU J073709. 1+653544 
detected in the 2004 August-October Chandra and XMM-Newton observations of the nearby spiral 
galaxy NGC 2403. The X-ray source exhibits X-ray pulsations with a period P~18 s and a nearly 
sinusoidal pulse shape and pulsed fraction 46-70% during the first three observations. A detailed timing 
analysis reveals a rapid decrease of the pulsation period from 18.25 s on Aug. 9 to 17.93 s on Sep. 12 
and possibly 17.56 s on Oct. 3, 2004. The X-ray spectra of CXOU J073709. 1+653544 are hard and are 
well fitted with an absorbed simple power law of photon index T ~ 0.9 — 1.2 in the 0.3-7 keV energy 
band. The X-ray properties of the source and the absence of an optical/UV counterpart brighter 
than 20 th magnitude allow us to identify CXOU J073709. 1+653544 as accreting X-ray pulsar in NGC 
2403. The maximum unabsorbed luminosity of the source in the 0.3-7 keV range, Ex ~ 2.6 x 10 38 ergs 
s" 1 at 3.2 Mpc, is at least 260 times higher than its quiescent luminosity, and exceeds the isotropic 
Eddington limit for a 1.4M compact object accreting hydrogen-rich material. The corresponding 
estimated luminosity in the 0.3-100 keV energy range could be as high as ~ 1.2 x 10 39 ergs s" 1 , 
assuming the typical pulsar energy spectrum with high-energy cut-off at 10-20 keV. The rate of 
decrease of the pulsation period of the source (P ~ — 10~ 7 s s~ x ) is one of the fastest observed among 
accreting pulsars. The evolution of the pulsation period suggests that it is dominated by the intrinsic 
spin-up of the compact object, which is almost certainly a neutron star. The X-ray luminosity of 
CXOU J073709. 1+653544 is high enough to account for the observed spin-up rate, assuming that 
the X-ray source is powered by disk accretion onto highly magnetized neutron star. Based on the 
transient behavior and overall X-ray properties of the source, we conclude that it could be an X-ray 
pulsar belonging to either a Be binary system or a low-mass system similar to the transient Galactic 
bursting pulsar GRO J1744-28. 

Subject headings: galaxies: individual (NGC 2403) — X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

The observations of accretion-powered X-ray pulsars 
provide a unique opportunity to study the dynami- 
cal properties of accreting matter and its interaction 
with the magnctosphcre of a neutron star, and the 
transfer of the angular momentum in a binary system 
(White, Swank & Holt 1983; Nagase 1989; Bildsten et 
al. 1997). Most of the known accreting X-ray pulsars 
are found in the high-mass binary systems involving a 
neutron star and OB donor star (either supergiant or 
Be type) and are associated with younger stellar popu- 
lations (Corbet 1986; Charles & Coe 2006). These two 
types correspond to different modes of accretion with 
supergiant systems accreting from a radially outflowing 
stellar wind, and the Be binaries accreting from a cir- 
cumstellar disk. Therefore, the accreting pulsars with 
supergiant companions are observed as persistent X-ray 
sources, while the Be systems are recurrent transients 
and highly variable sources and reach much higher X-ray 
luminosities up to 10 38 ~ 10 39 ergs s _1 . A small num- 
ber of X-ray pulsars is associated with low-mass binary 
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systems with mass donors ranging from a main-sequence 
stars to degenerate carbon-oxygen dwarfs and red giants 
and can reach persistent luminosities of > 10 38 ergs s _1 
(Nagase 1989; Bildsten et al. 1997). 

Until recently, the study of X-ray pulsars was limited to 
our Galaxy and the nearby Magellanic Clouds due to the 
limited sensitivity and spatial resolution of previous X- 
ray missions. With the advent of the Chandra and XMM- 
Newton X-ray observatories, it has become possible to 
study the spectral and timing properties of individual 
X-ray sources associated with more distant galaxies, in- 
cluding bright X-ray pulsators (Fabbiano & White 2006; 
Fabbiano 2006). For example, recent XMM-Newton ob- 
servations of M31 revealed a 865 s pulsating super- 
soft X-ray source and a 197 s accreting X-ray pul- 
sar candidate with luminosities of ~ 10 37 ergs s" 1 
and pulsed fractions of ~ 40% (Osborne et al. 2001; 
Trudolyubov et al. 2005). These results demonstrate 
that bright pulsating X-ray sources (Lx > 10 38 ergs s" 1 ) 
can be detected with XMM-Newton and Chandra in the 
galaxies beyond the Local Group up to the distances of 
a few Mpc. 

The spiral Scd galaxy NGC 2403 at 3.2 Mpc (Freedman 
& Madore 1988), provides a good opportunity to study 
X-ray source populations in a normal galaxy. NGC 2403 
was a target of previous X-ray observations with the Ein- 
stein (Fabbiano & Trinchieri 1987; Fabbiano et al. 1992), 
R OSA T (Roberts & Warwick 2000; Liu & Bregman 2005) 
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and ASCA observatories (Kotoku et al. 2000), which de- 
tected five discrete sources with one bright source clas- 
sified as an ultraluminous X-ray source (ULX) with lu- 
minosity above 10 39 ergs s _1 . Recent Chandra obser- 
vations covering the inner ~ 25% of the D25 area of 
NGC 2403 have revealed 41 discrete sources with appar- 
ent luminosities down to ~ 10 36 ergs s" 1 and allowed to 
study unresolved soft X-ray emission from the galactic 
disk (Schlegcl & Pannuti 2003; Fraternali ct al. 2002). 
In addition, the survey of luminous X-ray source popula- 
tions with XMM-Newton allowed to study spectral prop- 
erties of four brightest sources detected in NGC 2304 
(Winter et al. 2006) 

In this paper, we report on the discovery of the tran- 
sient pulsating X-ray source CXOU J073709. 1+653544 in 
NGC 2403 using publicly available August-October 2004 
observations made with Chandra and XMM-Newton as a 
part of the follow-up program studying recent supernova 
SN 2004dj (Pooley & Lewin 2004). 

2. OBSERVATIONS AND DATA REDUCTION 

In our analysis we used the data of four 2004 August- 
December Chandra/ ACIS-S observations of SN 2004dj 
region and 2001 April 17 observation of NGC 2403 (Ta- 
ble 1). The data of Chandra observations was processed 
using the CIAO v3.3 5 threads. We performed standard 
screening of the Chandra data to exclude time intervals 
with high background levels, and found no flares exceed- 
ing 20% of the average background level. For each ob- 
servation, we generated X-ray images in the 0.3-7 keV 
energy band, and used CIAO wavelet detection routine 
wavdetect to detect point sources. Only data in the 
0.3-7 keV energy range was used in the spectral anal- 
ysis of ACIS-S observations. The standard CIAO axbary 
tool was used to produce a barycenter-corrected source 
event lists used in the timing analysis. During three 2004 
August-October observations the source was positioned 
on the back-illuminated chip S3, offset by 1.3' ~ 1.4' from 
the aim point. To produce source filtered event lists and 
spectra, we extracted counts within the inclined ellipti- 
cal region with semi-axes of 2" and 2.5" centered at the 
position of CXOU J073709. 1+653544, large enough to 
include more than 90% of the source energy flux. Back- 
ground counts were extracted from the adjacent source- 
free regions with subsequent normalization by ratio of 
detector areas. 

We also use the data of 2003 April 30, 2003 Septem- 
ber 11 and 2004 September 12 XMM-Newton observa- 
tions of NGC 2403 (Table 1) with three European Photon 
Imaging Camera (EPIC) instruments (MOS1, MOS2 and 
pn)(Turner et al. 2001; Strueder et al. 2001), and the 
Optical Monitor (OM) telescope (Mason et al. 2001). We 
reduced XMM data using XMM-Newton Science Analysis 
System (SAS v 6.5.0) 6 . For 2004 September 12 observa- 
tion, we filtered the original EPIC event files based on 
the good time intervals produced from the light curves 
extracted from the whole detector in the 0.2-12 keV en- 
ergy band by applying an upper count rate threshold of 
6 and 2 counts s _1 for the EPIC pn and MOS cameras. 
The last ~ 20 ks of the EPIC observation is affected 
by high background, and has been excluded from anal- 

5 http://asc.harvard.edu/ciao/ 

6 See http://xmm.vilspa.esa.es/user 



ysis. The standard SAS tool barycen was used to per- 
form barycentric correction on the original EPIC event 
files used for timing analysis. A significant part of 2003 
April 30 and September 11 XMM-Newton observations 
was affected by background flares. After screening for 
high background only 4.7 ks (2003 Apr. 30 observation) 
and - 9 ks (2003 Sep. 11 observation) of EPIC-MOS 
exposure and 5 ks of EPIC-pn exposure (2003 Sep. 11 
observation) was left for the analysis (Table 1). 

We generated EPIC-pn and MOS images of NGC 
2403 field in the 0.3-7.0 keV energy band, and used 
the SAS standard maximum likelihood (ML) source 
detection script edetect_chain to detect point sources. 
We used bright X-ray sources with known counterparts 
from USNO-B catalog (Monet et al. 2003) and Chandra 
source lists to correct EPIC image astrometry. The astro- 
metric correction was also applied to the OM images, us- 
ing cross-correlation with USNO-B catalog. After correc- 
tion, we estimate residual systematic error in the source 
positions to be of the order 0.5 — 1" for both EPIC and 
OM. 

To generate EPIC-MOS source lightcurves and spectra 
during the 2004 September 12 XMM observation, we used 
a circular region of 18" radius centered at the position 
of CXOU J073709. 1+653544. Due to the source prox- 
imity to the edge of EPIC-pn CCD, the source counts 
were extracted from the elliptical region with semi-axes 
of 18" and 13" and position angle of 20 degrees, including 
~ 70% of the source energy flux. The adjacent source- 
free regions were used to extract background spectra and 
lightcurves. The source and background spectra were 
then rcnormalized by ratio of the detector areas. For 
spectral analysis, we used data in the 0.3 — 10 keV energy 
band. To synchronize source and background lightcurves 
from individual EPIC detectors, we used the identical 
time filtering criteria based on Mission Relative Time 
(MRT) , following the procedure described in Barnard et 
al. (2007). The background lightcurves were not sub- 
tracted from the source lightcurves, but were used later 
to estimate the background contribution in the calcula- 
tion of the source pulsed fractions. 

To estimate upper limits on the quiescent source lu- 
minosities, the ACIS and EPIC count rates were con- 
verted into energy fluxes in the 0.3-7 keV energy range 
using Web PIMMS 7 , assuming an absorbed power law 
spectral shape with photon index r = 1.0 and Galac- 
tic foreground absorbing column Nh=4 x 10 20 cm~ 2 
(Dickey & Lockman 1990). 

The energy spectra were grouped to contain a mini- 
mum of 20 counts per spectral bin in order to allow \ 2 
statistics, and fit to analytic models using the XSPEC 
v.12 8 fitting package (Arnaud 1996). EPIC-pn, MOS1 
and MOS2 data were fitted simultaneously, but with nor- 
malizations varying independently. For timing analysis 
we used standard XANADU/XRONOS v.5 9 tasks. 

In the following analysis we assume a distance of 3.2 
Mpc for NGC 2403 (Freedman & Madore 1988). All 
parameter errors quoted are 68% (lcr) confidence limits. 

3. RESULTS 

7 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 

8 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 

9 http://heasarc.gsfc.nasa.gov/docs/xanadu/xronos/xronos.html 
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3.1. Source Position 

We discovered a new X-ray source CXOU 
J073709. 1+653544 in the data of the 2004 August 
9 and detected it in the subsequent August 23 and 
October 3 Chandra observations of the NGC 2403 
field (Table 1). In addition, the analysis of archival 
2004 September 12 XMM-Newton observation revealed 
the presence of a bright X-ray source at the posi- 
tion consistent with CXO J073709. 13+653544. The 
source was not detected in the 2004 December 22 
Chandra observation (observation ID 4630). Com- 
bining the data of three Chandra observations, we 
measure the position of CXOU J073709. 1+653544 to be 
a = 07 ,l 37 m 09.139 s , S = +65°35'44.23" (J2000 equinox) 
with an uncertainty of ~ 0.5" (Fig. 1). The projected 
galactocentric distance of CXOU J073709. 1+653544 is 
- 2' or - 1.86 kpc at 3.2 Mpc. 

The search for the optical counterparts using the im- 
ages from the Digitized Sky Survey did not yield any 
stellar-like object brighter than bib ~ 20 within the 
error circle of CXOU J073709. 1+653544. We used 
the data of 2003 April 30 and 2004 Sept. 12 XMM- 
Newton/ 'OM observations to search for UV counterparts 
to the source prior and during its X-ray outburst (Fig. 
1). We did not detect any stellar counterparts to CXOU 
J073709. 1+653544 in the OM images down to the limit 
of - 20 m in the U and OM UVW1 (291 nm) bands. 

3.2. Long-Term Flux Evolution 

The evolution of X-ray luminosity of CXOU 
J073709. 1+653544 in the 0.3-7 keV energy band based 
on the data of 2004 Chandra and XMM-Newton obser- 
vations is shown in Fig. 2. The first two August 2004 
Chandra observations show the overall increase of the 
source X-ray luminosity in the 0.3-7 keV energy band 
from - 1.9 x 10 38 to - 2.3 x 10 38 ergs s" 1 . The subse- 
quent XMM-Newton and Chandra observations revealed 
the gradual decrease of the source luminosity to ~ 10 38 
ergs s _1 on Sept. 12, and ~ 5 x 10 37 ergs s" 1 on Oct. 3, 
2004. The source was not detected during the 2004 De- 
cember 22 Chandra observation with an upper limit on 
its 0.3-7 keV luminosity of < 10 36 ergs s" 1 (2a), > 230 
times lower than maximum measured outburst luminos- 
ity. It should be noted that due to the sparse sampling 
the lightcurvc of the source observed during 2004 August- 
October could be a result of either single long (> 60 days) 
outburst or the superposition of a number of shorter out- 
bursts poorly sampled in time. The analysis of archival 
2001 April 17 Chandra observation (observation ID 2014) 
(Schlegel & Pannuti 2003) of NGC 2403 taken prior to 
the 2004 outburst, allow us to estimate the average up- 
per limit (2a) on the source quiescent luminosity to be 
~ 10 36 ergs s _1 in the 0.3 — 7 keV energy band. The 
source was not detected in the two 2003 April 30 and 
September 11 XMM-Newton observations (Table 1), but 
with a much less stringent upper limit on its 0.3-7 keV 
luminosity, < 10 37 ergs s^ 1 (Fig. 2). 

3.3. X-ray Pulsations 

We performed a timing analysis of the CXOU 
J073709. 1+653544 using the data from C7iandra/ACIS-S 
and XMM-NewtonfEPIC detectors in the 0.3-7 keV en- 
ergy band. After a barycentric correction of the photon 



arrival times in the original event lists, we performed 
a Fast Fourier Transform (FFT) analysis using stan- 
dard XRONOS task powspec, in order to search for co- 
herent periodicities. For the analysis of XMM-Newton 
data, we used combined synchronized EPIC-pn and MOS 
lightcurves with 2.6 s time bins to improve sensitivity. 
We found strong peaks in the Fourier spectra of data 
from the 2004 Aug. 9 and 23 Chandra observations and 
2004 Sept. 12 XMM observation at the frequencies of 
(5.48-5. 58)xl0~ 2 Hz (Fig. 3, left panels). The strengths 
of the peaks in the individual Fourier spectra (Fig. 3) 
correspond to the period detection confidence of ~ 10 -10 , 
- lO" 11 and - 10~ 17 for the Aug. 9, 23 Chandra and 
Sept. 12 XMM- Newton observations (disregarding simul- 
taneous detection in all three observations). 

To estimate the pulsation periods more precisely, we 
used an epoch folding technique, assuming no period 
change during individual observations. The most likely 
values of the pulsation period (Table 2) were obtained fit- 
ting the peaks in the x 2 versus trial period distribution 
with a Gaussian. The period errors in Table 2 were com- 
puted following the procedure described in Leahy (1987). 
Then the source lightcurves were folded using the pe- 
riods determined from epoch folding analysis. The re- 
sulting folded lightcurves of CXOU J073709. 1+653544 in 
the 0.3-7 keV energy band during first three observations 
are shown in Fig. 3 (right panels). The source demon- 
strates quasi-sinusoidal pulse profiles in the 0.3-7 keV 
energy band during the first three observations (Fig. 3). 
The pulsed fraction, defined as (ImaxT m i n )/(I max +I min ), 
where I max and I m ; n represent source intensities at the 
maximum and minimum of the pulse profile excluding 
background photons, is significantly higher during the 
2004 Aug. 9 Chandra observation (70+4%) when com- 
pared to the observations #5 (46+4%) and #6 (56+3%) 
(Table 2). 

Although the analysis of the Oct. 3 Chandra obser- 
vation (obs. #7) did not show the presence of a strong 
single peak in the power density spectrum, we still per- 
formed a search for the pulsations in the data of this ob- 
servation using the epoch folding technique. The epoch 
folding analysis of the source lightcurve have led to the 
marginal detection of the pulsation with a period of 
17.565 s (Fig. 3, left panel). The folded lightcurve of the 
source in the 0.3-7 keV energy band with a significance 
of ~ 3.7a determined using % 2 for a nonvarying, constant 
model, is shown in the lower right panel of Fig. 3. The 
pulse profile possibly shows a non-sinusoidal structure, 
which probably explains why we do not see strong first 
harmonic in its Fourier spectrum. 

To investigate the energy dependence and time evolu- 
tion of the source pulse profile, light curves in the soft 
(0.3-2 keV) and hard (2-7 keV) bands were created for ob- 
servations #4, 5, and 6 10 and folded at the correspond- 
ing best pulsation periods (Figure 4). The correspond- 
ing background-corrected pulsed fractions in the soft and 
hard energy bands are given in Table 2. As can be seen 
from Fig. 4 and Table 2, the energy dependence of the 
pulse profile in the Aug. 9 Chandra observation (obs. 
#4) is qualitatively different from subsequent Chandra 

10 The limited source count statistics in the Oct. 3 observation 
does not permit a detailed study of energy dependence of pulsa- 
tions. 
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and XMM-Newton observations (obs. #5 and 6). The 
source modulation amplitude in the soft band is higher 
than in the hard band during Aug. 9 Chandra observa- 
tion, while the modulation is stronger in the hard band 
during the Aug. 23 Chandra and Sept. 12 XMM-Newton 
observations. The modulation fraction in the soft band 
is highest (78±5%) during the Aug. 9 observation, and 
drops to 43±6% and 50±6% on Aug. 23 and Sept. 12 
(Table 2). At the same time, the modulation fraction 
in the hard energy band remains compatible during all 
three observations (Table 2; Fig. 4). 

3.4. Long-term Period Evolution 

There is a significant change of the pulsation period 
of CXOU J073709. 1+653544 measured over the course 
of four observations covering a ~55-day period from 
P=18.253 s on Aug. 9, 2004 to P=17.934 s on Sept. 
12, 2004 and possibly P=17.565 s on Oct. 3, 2004 (Ta- 
ble 2; Fig. 5). During the first three observations cor- 
responding to a 35-day time interval, the source shows 
a roughly linear decrease of the pulsation period. Us- 
ing these three period measurements and assuming its 
linear time dependence, we estimate the rate of period 
change P ~ —(1.1 ± 0.2) x 10~ 7 s s -1 (pulse frequency 
derivative v ~ 3.4 x 10~ 10 Hz s -1 ), that corresponds 

to (P/P) 6 x 10~ 9 s -1 0.19 yr -1 for a period 

P = 18 s. 

3.5. X-ray Spectra 

The pulse phase averaged Chandra/ ACIS-S and XMM- 
Newton/EPlC spectra of CXOU J073709. 1+653544 can 
be adequately fit with the absorbed simple power law 
models with photon index, T ~ 0.9 — 1.2 and an equiv- 
alent hydrogen density Nh ~ (15 — 58) x 10 20 cm -2 . 
The corresponding absorbed luminosity of the source 
in the 0.3-7 keV band ranges between ~ 5 x 10 37 and 
~ 2.3 x 10 38 ergs s _1 , assuming the distance of 3.2 Mpc. 
The best-fit spectral model parameters for individual ob- 
servations are given in Table 2. The luminosity and en- 
ergy spectra of the source are similar to that of the bright 
X-ray pulsars in our Galaxy and the Magellanic Clouds 
(Nagase 1989; Yokogawa et al. 2003). The shape of the 
energy spectra of CXOU J073709. 1+653544 in the 0.3-7 
keV energy band remains essentially constant throught 
the outburst despite significant change of the source lu- 
minosity. For all observations, the measured absorbing 
column Nh is ~3-10 times higher than the Galactic hy- 
drogen column in the direction of NGC 2403, 4xl0 20 
cm -2 (Dickey & Lockman 1990), consistent with an ad- 
ditional intrinsic absorption within the system and inside 
the disk of NGC 2403. 

4. DISCUSSION 

A combination of transient behavior and the extremely 
large spin-up rate of CXOU J073709. 1+653544 makes it 
a unique object among pulsating X-ray sources. The 
overall X-ray properties (X-ray spectrum, pulsation pe- 
riod etc.) of the source are consistent with that of the 
accreting X-ray pulsar in a binary system (White, Swank 
& Holt 1983; Nagase 1989). The absence of bright opti- 
cal counterparts (both during quiescence and outburst), 
transient behavior, overall X-ray properties of CXOU 
J073709.1+653544 and positional coincidence with NGC 



2403, allow us to conclude that is almost certainly located 
outside our Galaxy and belongs to NGC 2403. Placing 
CXOU J073709. 1+653544 at the distance of NGC 2403 
(3.2 Mpc) implies very high luminosities of the source 
during the outburst ((0.6 ~ 2.4) x 10 38 ergs s -1 in the 
0.3-7 keV energy band and (0.8 ~ 3.8) x 10 38 ergs s -1 
in the 0.3-10 keV energy band) with maximum observed 
luminosity exceeding the isotropic Eddington luminosity 
limit for a 1.4M Q object accreting hydrogen- rich mate- 
rial. If CXOU J073709. 1+653544 is indeed an accreting 
X-ray pulsar in NGC 2403, the relative rate of change 
of its pulsation period, (P/P) ~ —6 x 10~ 9 s" 1 is the 
highest observed in this class of objects to date (Nagase 
1989; Bildsten et al. 1997). 

If the X-ray pulsations observed in CXOU 
J073709. 1+653544 result from rotation of a highly 
magnetized accreting compact object in a binary sys- 
tem, the change of its pulsation period can be explained 
by combination of binary orbital motion of the pulsar 
(Skinner et al. 1982) and transfer of angular momentum 
in the accretion process (Rappaport & Joss 1977; Ghosh 
& Lamb 1979). The relatively long time scale on which 
we observe the decrease of the pulsation period and 
its high, nearly constant rate suggest that the intrinsic 
spin-up of the pulsar should make the major contri- 
bution to the observed change of the pulsation period. 
Using the observed pulsation parameters of CXOU 
J073709. 1+653544, and assuming that disk accretion is 
occurring in the system, one can make a simple estimate 
of the lower limit on the accretion rate required to 
explain the observed spin-up rate P. Let us assume 
that all the specific angular momentum of the accreting 
material is transferred to the compact object at the mag- 
netospheric radius r m . In order to keep the "propeller" 
effect from preventing accretion onto compact object, 
the magnetospheric radius must be less that the corota- 
tion radius r m < r c = (GM x /4tt 2 ) 1/3 P 2/3 , where P is 
the rotational period and Mx is a compact object mass 
(Illarionov & Sunyaev 1975). The maximum torque 
experienced by compact object N — M(GMxr c ) 1/>2 , 
where M denotes the accretion rate through the disk, 
should be greater or equal to the observed torque 
N bs = Id) = 2wI(P/ P 2 ), where / is the moment of 
inertia of the compact object. As a result, we obtain 

M >I(47r 2 /GMx) 2/3 PP~ 7/3 (1) 

or, 

M>5.2xlO'v( TjF ^ r ) (£f'V\ (2) 

where / = (J/10 45 gem 2 ) (M X /M Q )~ 2 / 3 . Substitut- 
ing typical white dwarf parameters (/ ~ 10 50 g cm 2 , 
Mx = Mq) and the observed rotation period and its 
derivative into equation (2) gives an unreasonably high 
mass accretion rate M > 5 x 10 23 g s -1 . On the other 
hand, for a typical neutron star parameters (/ ~ 10 45 g 
cm 2 , Mx = 1.4M©) the estimated mass accretion rate 
required to produce the observed spin-up is M > 4 x 10 18 
g s -1 = 6.3 x 10~ 8 Mq year -1 . The corresponding bolo- 
metric luminosity of the IAMq neutron star accreting at 
such rate should be Lboi > 7.3 x 10 38 ergs s _1 , assuming 
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that all gravitational energy of the infalling matter re- 
leased in the process of accretion is converted into radia- 
tion. The estimated bolometric luminosity of the source 
is very high, exceeding the Eddington critical luminos- 
ity for a 1.4M Q star accreting hydrogen-rich material at 
least by factor of 4. 

It is interesting to compare the luminosity required 
from accretion spin-up mechanism to the total X-ray 
luminosity of the source during the outburst. In or- 
der to estimate the total X-ray luminosity of CXOU 
J073709. 1+653544, we extrapolated the results of our 
spectral analysis to the higher energies, assuming that 
the broad-band spectrum of the source in the 0.3-100 kcV 
band is represented by a power law with high-energy cut- 
off at energies above 10-20 keV, as typically observed in 
the accreting X-ray pulsars (White, Swank & Holt 1983; 
Nagase 1989). To model the broad-band spectrum, we 
used a simple power law model with photon index and 
normalization derived from spectral fitting in the 0.3-7 
keV band, modified at energies above a high-energy cut- 
off E c by the function exp[(E c — E)/Ef], where Ef is the 
folding energy. For the Aug. 23 observation, the esti- 
mated X-ray luminosity of CXOU J073709. 1+653544 in 
the 0.3-100 keV energy band falls between - 8 x 10 38 
and 1.2 x 10 39 ergs s" 1 for E c = 10 ~ 20 keV and 
Ef = 10 keV. Note that this estimate should be regarded 
as a lower limit, as it does not account for the possible 
presence of the soft (kT ~ 0.1 — 0.2 kcV) component in 
the spectrum of the source, similar to that observed in 
some of the bright X-ray pulsars (Yokogawa et al. 2003). 
Therefore, the estimated maximum X-ray luminosity of 
CXOU J073709. 1+653544 is high enough to account 
for the observed spin-up rate, assuming that the X-ray 
source is powered by disk accretion onto highly magne- 
tized neutron star. 

The estimated maximum luminosity of the source can 
be used to put an upper limit on the strength Bq of the 
magnetic field of the neutron star. Assuming the dipolc 
configuration of the pulsar magnetic field, and requiring 
the magnetospheric radius to be smaller than corotation 
radius r m < r c , we obtain 

where h = (R/10 6 cm)- 5 / 2 (M^/M©) 1 / 3 and R is the 
radius of the neutron star. 

Although both the observed luminosity of the source 
and luminosities required by spin-up due to disk accre- 
tion are extremely high, they lie within the range ob- 
served for pulsating neutron stars that can reach lumi- 
nosities of <~ 5 — 10 x 10 38 ergs s _1 (White & Carpenter 
1978; Skinner et al. 1982; Nagase 1989). The required lu- 
minosities are also consistent with theoretical predictions 
for super-Eddington accretion onto highly magnetized 
(B > 10 12 G) neutron star (Basko & Sunyaev 1976). At 
high accretion rates, the accreting matter is expected to 
release its gravitational energy inside accretion columns 
above the magnetic poles of the neutron star. Since most 
of radiation is emitted from the sides of these columns, 
it can escape the system without crossing the bulk of 
the inflowing accreting matter. This anisotropy allows 
the luminosity of the neutron star to exceed formal Ed- 
dington limit by large factors. For the accretion column 



with a cross-section in the form of a thin ring segment 
of length I and thickness d, the limiting luminosity that 
can be reached in the Basko & Sunyaev (1976) model is 
given by 

«■»(*£) (4) 

Therefore, for the highly super-Eddington accretion 
rate and I ^> d the emission from accretion columns 
could explain the observed high luminosity of CXOU 
J073709. 1+653544. 

The measured pulse period places CXOU 
J073709. 1+653544 among Be systems or short-period 
binary supcrgiant systems with Roche lobe overflow 
on a Corbet diagram (pulse period-orbital period 
plot) (Corbet 1986; Waters & van Kerkwijk 1989; Bild- 
sten et al. 1997). The transient behavior of the source 
suggests a higher probability for it to be a member of 
Be binary class, since majority of Be systems display 
recurrent /transient outbursts. The high peak luminosity 
of CXOU J073709. 1+653544 (L x ~ 10 39 ergs s" 1 ) is 
also not unusual for the Be X-ray pulsars with transient 
Be X-ray pulsar A0538-66 reaching similar luminosities 
during its giant outbursts (Skinner et al. 1982). 

Using the Corbet diagram, one can try to predict the 
value of the binary period of CXOU J073709. 1+653544. 
If it is indeed a Be system, the estimated orbital period 
P or b can vary between ~ 20 and ~ 100 days. The ob- 
served duration of the source outburst allows us to put 
further constraints on the value of the binary period. 
Assuming that we observed a single outburst from the 
source and that it corresponds to one orbital cycle, wc 
obtain 60 < P or b < 100 days. Note, however, that if the 
source went through the giant (class II) outburst lasting 
for several orbits (Stella, White & Rosner 1986) or a se- 
quence of shorter outbursts during August-October 2004, 
the orbital period of the system could be much shorter 
than 60 days. 

Another possibility is that CXOU J073709. 1+653544 
is a transient low-mass binary system. The transient 
behavior and high luminosity of the source are some- 
what similar to that of the Galactic bursting X-ray pulsar 
GRO J1744-28 (Kouveliotou et al. 1996), that reached a 
luminosity of ~ 10 39 ergs s _1 during the peak of its out- 
burst (Giles et al. 1996; Sazonov, Sunyaev & Lund 1997). 
In the low-mass system scenario, the transient outbursts 
of the source could be explained by sudden increase of 
accretion from the disk reservoir (Sunyaev & Shakura 
1977) or as a result of the viscous-thermal instability of 
the quiescent accretion disk. 

The discovery of X-ray pulsations in CXOU 
J073709. 1+653544 is yet another demonstration 
that deep monitoring observations with Chandra and 
XMM-Newton have a great potential in revealing 
populations of both persistent and transient pulsating 
X-ray sources in nearby galaxies. Future observations of 
CXOU J073709. 1+653544, if it reappears, could possibly 
give deeper insight into the behavior of this and other 
similar X-ray pulsar systems. The complete coverage of 
the source outbursts with more frequent sampling would 
allow to estimate the orbital parameters of the system 
and help to determine its nature. 
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TABLE 1 

Chandra and XMM-Newton Observations of NGC 2403 Used in the Analysis. 



Obs.# 


Date, UT 


Date, T.TD 


Obs. ID 


Mission/Instrument 


Mode/ 


RA (J2000) a 


Dec (J2000) a 


Exp> 














Filter 


(h:m:s) 


(d:m:s) 


(ks) 


1 


2001 Apr. 


17 


12016.401(0.221) 


2014 


CTiandra/ ACIS-S 


Faint 


07:36:49.86 


+ 65:34:43.5 


35.5 


2 


2003 Apr. 


30 


12759.756(0.090) 


0150651101 


XMM-Newtm/EPIC 


Full/Thin 


07:36:51.00 


+ 65:36:10.0 


4.7(MOS) c 


3 


2003 Sep. 


11 


12893.044(0.065) 


0150651201 


XMM-Newton/EPIC 


Full/Thin 


07:36:51.00 


+ 65:36:10.0 


8.9(MOS)/5.0(pn) 


4 


2004 Aug. 


09 


13226.561(0.273) 


4627 


Chandra/ ACIS-S 


Faint 


07:37:22.34 


+ 65:36:10.7 


40.0 


5 


2004 Aug. 


23 


13240.701(0.290) 


4628 


Chandra/ ACIS-S 


Faint 


07:37:21.72 


+ 65:36:10.9 


46.5 


6 


2004 Sep. 


12 


13261.182(0.449) 


0164560901 


XMM-Newton/EPIC 


Full/Modiur 


a 07:37:17.00 


+ 65:35:57.8 


59.0(MOS)/52.4(pn) 


7 


2004 Oct. 


03 


13281.951(0.276) 


4629 


Chandra/ ACIS-S 


Faint 


07:37:19.45 


+65:36:29.8 


44.5 


8 


2004 Dec. 


22 


13361.938(0.304) 


4630 


Chandra/ ACIS-S 


Faint 


07:37:12.85 


+ 65:36:20.9 


49.9 


a pointing coordinates 



"instrument exposure used in the analysis 

c the source position coincides with the gap between EPIC-pn CCDs, only EPIC-MOS data was used in the analysis 
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TABLE 2 

X-ray Pulsation Parameters and Power Law Spectral Fit Information for CXOU J073709. 1+653544. 



Obs. 


Period 

« 


PF 0.3-7fceV 

(%)» 


PF 0.3-2feeV 

(%) a 


PF 2-7fce 
(%) a 


V N H 

(X10 20 cm^ 2 ) 




Photon 
Index 


Flux b 




Flux c 




x 2 

(d.O.f) 






Instrument 


4 


18.25247(73) 


70 ± 4 


78 ± 5 


63 ± 6 


37+« 


1 


18 -0.12 


1.58 ± 0.07 


1 


96 ± 0. 


12 


25.6(25) 


1.94 


2.40 


ACIS-S 


5 


18.09825(53) 


46 ± 4 


43 ± 6 


56 ± 5 


20 ± 5 


0. 


96 ± 0.10 


1.90 ± 0.07 


2 


12 ± 0. 


14 


27.7(33) 


2.33 


2.60 


ACIS-S 


6 


17.93443(30) 


51 ± 3 


50 ± 6 


64 ± 4 


«t; 


0. 


88 ± 0.08 


0.89 ± 0.04 





99 ± 0. 


06 


52.4(64) 


1.09 


1.21 


EPIC/pn + MOS 


7 


17.5645(14) 








»±3 


1 


■^-0.34 


0.38 ± 0.03 





50 ± 0. 


06 


3.0(4) 


0.47 


0.61 


ACIS-S 



7, 0.3 



r gy 



ichncd as 



source background-corrected intensities at the maximum and minimum of the pulse profile 



+lmin)> where i max and l min represent 
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b absorbed model flux in the 0.3 — 7 keV energy range in units of 10 
c unabsorbed model flux in the 0.3 
d absorbed luminosity in the 0.3 

e unabsorbed luminosity in the 0.3 — 7 keV energy range in units of 10 38 



erg s 



cm 

-l 



7 keV energy range in units of 10 erg s " < i t 
7 keV energy range in units of 10 erg s 



erg s 



assuming the distance of 3.2 Mpc 
, assuming the distance of 3.2 Mpc 



Chandra/ACIS-S; Aug. 23 2004; 0.3-7.0 keV 



XMM-Newton/EPIC-MOS; Sept. 12 2004; 0.3-7.0 keV 



SN 2004dj 



CXOU J073709.1 +653544 



S 36:00.0 



% 35:00.0 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 



CXOU J073709.1 +653544 



7 

SN 2004dj 




24.0 12.0 7:37:00.0 

R.A. (J2000) 

XMM-Newton/OM; Sept. 12 2004; UVW1 



_i_ 



i .... i 



_i_ 



48.0 36.0 24.0 12.0 7:37:00.0 48.0 36.0 36:24.0 
R.A. (J2000) 

NGC 3403; DSS-2 



SN 2004dj 



CXOU J073709.1 +653544 




J » 



10" 




18.0 16.0 14.0 12.0 10.0 08.0 06.0 04.0 02.0 7:37:00.0 
R.A. (J2000) 



3D. II 



7J7:DD.Q 3D.D 
R.A. :'J£000.. 



Fig. 1. — Upper left: Chandra/ AC1S-S image covering a 5' X 5' region of NGC 2403 centered on the CXOU J073709. 1+653544, taken on 
Aug. 23, 2004. Upper right: Combined 0.3-7 keV XMM/EPIC-MOS image covering a 10' X 10' region of NGC 2403 centered on the CXOU 
J073709. 1+653544 position, taken on Sept. 12, 2004. In both upper panels the positions of a new pulsar CXOU J073709. 1+653544 and 
supernova SN 2004dj are marked with arrows, and the position of NGC 2403 nucleus is shown with a cross. Lower left: XMM-Newton/OM 
UVW1 band image of NGC 2403 field taken on Sept. 12, 2004. The image is a 2' X 2' square centered on the CXOU J073709.1+653544 
position. The localization of CXOU J073709. 1+653544 is shown with black circle of 1.5" radius (3<r). Lower right: Optical B-band image 
of NGC 2403 from the Second Generation Digitized Sky Survey. The position of CXOU J073709.1+653544 is shown with white circle. The 
regions covered by Chandra/ ACIS, XMM/EPIC and XMM/OM images are shown with black, gray and white squares. 
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CXOU J073709. 1+653544 
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Fig. 2.— (Upper panel:) The X-ray light curve of CXOU J073709. 1+653544 in the 0.3-7 keV energy band based on the results of 
2004 Chandra and XMM-Newton observations. (Lower panel:) long-term light curve of CXOU J073709. 1+653544 based on the results of 
2001-2005 Chandra and XMM-Newton observations. The luminosities were calculated assuming a source distance of 3.2 Mpc. 
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Fig. 3. — (Left upper three panels) Power spectra of CXOU J073709. 1+653544 obtained using the data of 2004 Aug. 9 (upper panel) and 
Aug. 24 (middle panel) Chandra/ ACIS-S observations and 2004 Sept. 12 XMM-Newton/EPIC observations (lower panel) in the 0.3-7 keV 
energy band. (Right upper three panels) Corresponding pulse profiles folded with most likely pulsation periods. The background levels are 
represented by the dotted lines. The ephemeris is defined arbitrarily such as the dip in the pulse profile falls at phase 0.8. (Lower panels) 
The results of epoch-folding analysis of the data of 2004 Oct. 3 Chandra/ ACIS observation of the source. 
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Aug ; 9 3004; Chandra/ACIS-S Aup. 33 3004; Chandra/ACIS-S Sept 13, 3004; XMM/EPIC 




0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 

Pulse Phase Pulse Phase Pulse Phase 

Fig. 4.— X-ray lightcurvcs of CXOU J073709. 1+653544 during the 2004 Aug. 9 and 23 Chandra and Sept. 12 XMM-Newton observations 
folded at the corresponding best periods (Table 2) in the soft (0.3-2 keV) and hard (2-7 keV) energy bands (upper middle panels) along 
with hardness ratios, computed taking background contribution into account (bottom panels). The background levels are represented by 
the dotted lines. 
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Fig. 5. — The long-term evolution of the pulsation period of CXOU J073709. 1+653544. No orbital corrections have been applied to the 
pulse frequencies, since the orbital parameters are unknown. 
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Fig. 6. — Left panels: Count spectra and model ratios of CXOU J073709. 1+653544 obtained during August 23 Chandra/ ACIS-S 
observation. Right panels: EPIC-pn count spectra and model ratios of the source during September 12 XMM-Newton observation. The 
best-fit absorbed power law model approximation is shown with thick histograms. 



